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Abstract 

Neonicotinoid insecticides act on nicotinic acetylcholine receptor and are particularly effective against sucking pests. They 
are widely used in crops protection to fight against aphids, which cause severe dannage. In the present study we evaluated 
the susceptibility of the pea aphid Acyrthosiphon pisum to the commonly used neonicotinoid insecticides imidacloprid (IMI), 
thiamethoxam (TMX) and clothianidin (CLT). Binding studies on aphid membrane preparations revealed the existence of 
high and low-affinity binding sites for pH]-IMI (Kd of 0.16±0.04 nM and 41.7±5.9 nM) and for the nicotinic antagonist [^^^1]- 
a-bungarotoxin (Kd of 0.008±0.002 nM and 1.135±0.213 nM). Competitive binding experiments demonstrated that TMX 
displayed a higher affinity than IMI for [^^^l]-o(-bungarotoxin binding sites while CLT affinity was similar for both [^^^l]-oi- 
bungarotoxin and pH]-IMI binding sites. Interestingly, toxicological studies revealed that at 48 h, IMI {LC5o = 0.038 |xg/ml) 
and TMX (LC5o = 0.034 |xg/ml) were more toxic than CLT (LC5o = 0.118 ng/ml). The effect of TMX could be associated to its 
metabolite CLT as demonstrated by HPLC/MS analysis. In addition, we found that aphid larvae treated either with IMI, TMX 
or CLT showed a strong variation of nAChR subunit expression. Using semi-quantitative PCR experiments, we detected for 
all insecticides an increase of ApisumalO and Apisumpi expressions levels, whereas ApisumP2 expression decreased. 
Moreover, some other receptor subunits seemed to be differently regulated according to the insecticide used. Finally, we 
also demonstrated that nAChR subunit expression differed during pea aphid development. Altogether these results 
highlight species specificity that should be taken into account in pest management strategies. 
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introduction 

Neonicotinoid insecticides include several compounds such as 
imidacloprid (IMI), clothanidin (CLT) and thiametoxam (TMX). 
They are efficient agonists of insect neuronal nicotinic acetylcho- 
line receptors (nAChRs) which are pentameric receptors formed 
by identical (homomeric) or diflFerent (heteromeric) subunits 
[1,2,3]. In the context of a stronger legislation on insecticide use 
to limit environmental and health concerns, approaches to 
describe and understand the cellular and molecular mechanisms 
involved in insecticide resistance are needed. Insects represent a 
very diverse group of animals and most Orders diverged 
approximately 300 million years ago [4]. As a consequence, 
adaptive mechanisms that confer insecticide resistance can vary 
from one order/species to another. Among insect pests that cause 
damage to agriculture, aphids (Hemiptera) have a particular 
biology: they feed from phloem sap (not by chewing plant tissues). 



and thus transmit plant viruses very efficiendy [5]. Their pest 
status is also attributable to their pecuhar reproductive mode [5] . 
Asexual reproduction of aphids by parthenogenesis (during spring 
and summer) leads to extremely rapid population growth [6]. 
Several insecticides, such as neonicotinoids, are used as seed 
treatment to limit the impact of sucking-pest like aphids, Sitohion 
avenae, Aphis cmccivora and Myzus perdcae [7,8]. 

Bioassay studies have revealed that neonicotinicotinoid suscep- 
tibility varies between insect species [8,9,10]. In Hemiptera, IMI 
showed susceptibility difTerences with LC95 values between 0.32 
and 40 mg.L~' [9]. Moreover, studies performed with the cotton 
aphid Aphis gossipii demonstrated that one IMI resistant-strain was 
still susceptible toward TMX and CLT suggesting no cross 
resistance [10] whereas a cross resistance was found in AI)izus 
persicae, with resistance factors of 1 1 , 18 and 1 00 for IMI, TMX 
and CLT, respectively [8]. This discrepancy suggested that some 
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aphid species could carry particular resistance mechanisms. 
Similar difTerences could be found using competitive binding 
studies. In Aphis craccivora it was demonstrated that TMX was a 
non-competitive inhibitor for [^H]-IMI suggesting that it binds to 
a different site or in a different mode than IMI while in Myzfis 
persicae, it was found that the resistant strain with the R81T 
mutation on the MpPl subunit developed cross-resistance against 
IMI and TMX demonstrating that they interact with the same site 
[3]. In addition, saturation studies suggest that IMI binds to high- 
and low-aflinity binding sites in the aphids M. persicae and A. 
craccivora [11,12]. Similar high and low affinity binding sites were 
also itlc'ntificd in M. persicae for the nicotinic antagonist, ot- 
Bungarotoxin (a-Bgt) [13]. Altogether these data demonstrati; that 
aphid species can carry different sensitivities against neonicotinoid 
insecticides. 

In this study, we evaluated the binding properties and 
toxicological effects of IMI, TMX and CLT on the pea aphid 
and showed that IMI and TMX bind two different nAChR 
populations, both binding CLT, although the toxicity of CLT is 
lower than toxicities of IMI and TMX. In addition, using the fuU 
sequenced genome of A. pisum [14], we demonstrated that the 
eleven A. pisum nAChR subunits previously identified [15] are 
differentially regulated during aphid development and after 
neonicotinoid intoxication. 

Materials and Methods 

Insects 

The pea aphid (A.pisurri) sequenced strain LSRl (corresponding 
to the reference genome) was generously provided by INRA- 
Rennes IGEPP. Unwinged parthenogenetic females were reared 
on faba bean [Vicia fabae) plants in a 16L: 8D photoperiod at 
constant temperature of 22°C in a climate chamber. Under these 
conditions, aphids reproduce by viviparous parthenogenesis, as 
clonal female: new born larvae become adults after four molts. 
Adults and larvae at each of the 4 stages were collected. Larval 
stages were determined by identifying the number of antennal 
segments. 

Insecticides 

IMI, TMX and CLT were purchased from Sigma-Aldrich 
(Saint-Louis, USA). Insecticides were dissolved in DMSO to give 
final concentrations of 50 mg.ml"'. For binding experiments, 
insecticides were dissolved in the corresponding buffer (PBS buffer 
or Tris-HCl buffer for ['^'^II-a-Bgt and [^H]-IMI experiments 
respectively). For intoxication experiments, insecticides were 
dissolv(;d in artificial diet at a final concentration of 0.2% DMSO 
for lOOng/ml and 2% DMSO for 1000 ng/ml. These concentra- 
tions of DMSO were used as controls in insecticidal assays. 

Binding assays 

Membrane preparations were isolated from frozen aphids 
according to the Wiesner and Kayser protocol [12]. Whole aphids 
were homogenized with a pestle motor in 4°C dissociation 
medium at pH 7.0. The dissociation medium contained: 20 mM 
sodium phosphate, 150 mM sodium chloride, 1 mM EDTA, 
0.1 mM phenylmethanesulfonyl fluoride (dissolved in acetone), 
and 2 (tg each of pepstatin, chymostatin, and leupeptin (dissolved 
in methanol, DMSO, and water, respectively). After homogeni- 
zation, samples were centrifuged 10 min at 1000 g and superna- 
tant was collected and ultracentrifuged 30 min at 4°C and 
43000 g. The precipitate was washed with cold dissociation 
medium and then ultracentrifuged. The final pellet was resus- 
pended in 3 ml of dissociation medium. Total protein was 



quantified by spectrofluorometry at 750 nm according to the 
Lowry colorimetric method (DC protein assay, Biorad, France) 
with a range of BSA as a standard. Membrane preparations were 
conserved at — 80°C until use. 

a-Bgt binding experiments were performed using 40 |j,g of 
aphid total membrane protein in a total volume of 300 |J,1 of PBS 
Buffer (Na2HP04, NaH2P04, NaCl, pH= 7.2) + 0.1% of bovine 
serum albumin and ['^'I]-0(-Bgt (2200 Ci/mmol, PerkinElmer, 
USA) as radiolabeled ligand for total binding measurement. For 
non-specific binding determination, 1 |j,M of cobratoxin was 
added prior to membrane incubation. In saturation assays, the 
concentration of ['^^I]-a-Bgt varied from 18 nM to 0.9 pM to 
obtain a complete saturation binding curve. Competitive assays 
were performed with IMI, TMX and CLT. For this purpose, 
membranes were incubated with various concentrations of 
unlabeled competitor and ['^^I|-a-Bgt at 0.08 nM for CLT and 
TMX and 0.6 nM for IMI, respectively. Incubations were 
performed at room temperature during 4 hours and terminated 
by rapid vacuum filtration using GF/C Glass microfiber filters 
presoaked in polyethyleneiminine 0.5%. Filters were rapidly 
washed (< 20 s) twice with 5 ml of cold PBS Buffer at 0.01 M 
and transferred in tubes for immediate counting on a y-counter. 

IMI binding experiments were performed using 200 Jig of aphid 
total membrane protein in a final volume of 300 \A of Tris-HCl 
Buffer (10 mM, pH = 7.4) and [^H]-IMI (ARC, 40 Ci/mmol) as 
radiotracer for total binding measurement. For non-specific 
binding determination, 0.3 mM of unlabeled IMI was added 
prior to membrane incubation. In saturation assays, the concen- 
tration of ['^H]-IMI varied from 500 nM to 0.5 pM to obtain a 
complete binding curve. Competitive assays were achieved for 
CLT, TMX, IMI and a-Bgt using 25 nM of ['^H]-IMI. 
Incubations were performed at room temperature during 4 hours 
and terminated by rapid vacuum filtration using GF/C Glass 
microfiber filters presoaked in 0.5% polyethyleneiminine. Filters 
were rapidly washed (< 20s) with cold Tris-HCl Buffer and dried 
for 1 hour before incubation in 5 ml of scintillation hquid 
(PerkinElmer, USA) and counting. 

Insecticidal assays 

The susceptibility of A. pisum to IMI, TMX and CLT was 
determined using an artificial diet bioassay according to Sadeghi et 
a/. [16]. In brief, adults were put on a feeding apparatus (day D-1) 
containing 200 |J,1 of artificial diet. Then, the first-instar nymphs 
were transferred (day DO) to freshly prepared diet with insecticide 
added (treatment series) or DMSO added (control series). For each 
insecticide, eight concentrations ranging from 0.001 to 
1000 |ig.ml were tested. The mortality was scored after 24 h 
(day D+1) and 48 h (day D+2). Aphids that were unable to walk 
were considered dead [8] and were removed. Corrected mortality 
percentages were calculated using Henderson Tilton's formula 
after 24 h (day D-l-1) and 48 h (day D-l-2) of insecticide exposure. 

HPLC-MS/MS analysis 

For tissue extraction, 2.4 g of TMX (at LC50) treated aphid 
larvae were extracted with an acidified (0.2% acetic acid) water- 
methanol mixture (50/50) in a 50 ml centrifuge tube. Sample 
extracts were then filtered and purified on Oasis HLB Cartridges 
(200 mg) (Waters SAS, France). Elution was performed using 6 ml 
of acetonitrUe. The obtained acetonitrUe extract after the elution 
of the Oasis cartridges was reduced to dryness and the residue was 
re-dissolved in 1 ml methanol-water mixture (10/90). 40 |J.l of 
each sample was analyzed by high-performance liquid chroma- 
tography coupled to tandem mass spectrometry (HPLC-MS/MS). 
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HPLC-MS/MS was performed with an Ultimate 3000 rapid 
separation liquid chromatography system (Dionex, USA) coupled 
to an API 4000 Qtrap MS/MS from Applied Biosystems (Foster 
City, CA, USA). Separation was performed on a Phenomenex 
(Torrance, CA, USA) C18 column at 35°C with a gradient of 
water/methanol/acetic acid at a flow rate of 0.2 ml.min and 
5 mM ammonium acetate. MS/MS detection was performed in 
the multi-reaction-monitoring (MRM) mode using an ESI 
interface in the positive ion mode. The ionization voltage was 
5500V, and the nebulizer and curtain gases were at 50 psi and 
25 psi, respectively. The drying gas to assist the solvent 
evaporation in the source (600°C) was at 40 psi. 

Optimisation of MRM transitions, collision energies and cone 
voltage were performed by direct injection of standard solutions. 
The optimized parameters for the detection of the two compounds 
(TMX and CLT) are listed in Table 1. With these parameters, 
calibration curves were linear over the concentration range of 0.9 
to 20 |J,g.l with a correlation coefficient (r) greater than 0.99. 

Expression of nAChR subunits during developmental 
stages 

Total RNAs were extracted from A. pisum adults, at different 
larval stages (Stage LI to L4), using RNA Easy mini Plant Kit 
(Qiagen, Courtaboeuf France). To avoid genomic DNA (gDNA) 
contamination, total RNAs were treated using DNAse I kit 
(Invitrogen, Carlsbad, USA) according to manufacturer recom- 
mendations. RNAs were retro-transcribed using random hexamers 
with RevertAid kit (Thermoscientific, Waltham, USA), dissolved 
in RNAse-free water and conserved at — 20°C. DNAse treatment 
was validated by PCR using primers set amplifying intron- 
containing sequence. Primer sets (Table 2) were designed using 
PrimerS software based on the A. pisum genome (http://www.ncbi. 
nlm.nih.gov/genbank/). Because Dale et al. identified potential 
alternatively spliced isoforms for Apisuma4 (exon 4) ApisumotG 
(exon 6) and Apisumo(7 (exon 6 and 7), primers were designed out 
of these exons, using genome information [14]. Amplification 
specificity of each primer set was also verified by cloning and 
sequencing the amplification products (data not shown). Amph- 
fication efficiencies were between 88 and 109%, allowing 
validation of each primer set for qPCR experiments. Because 
none of the endogenous reference genes had stable expression 
during developmental stages, external reference gene, luciferase, was 
used for normahzation, as previously described [17,18]. Thus, 
10 pg/ 1000 ng of luciferase RNA (Promega, Fitchburg, Wisconsin 
USA) were added after RNA extraction [17,18]. 

qPCR experiments were optimized according to MIQE 
Guideline recommendations [19] using ABI Prism 7700 instru- 
ment and 2X SYBR Green PCR Master Mix (Applied Biosystems, 
Courtaboeuf, France). Experiments were performed in triplicate 
using 100 ng of total RNA and 150 nM primers in a final volume 
of 25 |il. Product specificity was further assessed by electrophoresis 
on a 2% agarose gel with a 50 bp ladder and by dissociation 
curves giving rise to a single peak at the specific melting 
temperature [20]. Relative expression ratio (R) was calculated 
according to the Pfaffl formula [21], using primer efficiency (E) 
and CP value variation between control and sample (ACP) for 
each nAChR subunit. Ratio were normalized to reference genes 
and expressed in percentage: R = (E.^b,,,;,)"^^^™''"""'""*"' ' '™p''=> 

/ /T7 \ AC^Preference fcontrpl — sample) T 'X" 1 

/ (Ereference) ■ Lucuerasc was uscd as 

reference gene and quantification was relative to the first larval 
stage (LI). 
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Table 2. Primers used to amplify nicotinic acetylcholine receptors subunlts in quantitative PCR experiments. 



gene 


forwBrd primer 




reverse primer 




Size (bp) 




name 


nucleotidic sequence 


name 


nucleotidic sequence 




a 1 


qpAlSl 


CGGTCATTGTCGGTCAGTTG 


qpAIRl 


TGGCATCGGCACTTCCAT 


60 


a 2 


qpA2S2 


GGTCGTCACCATCATCATC 


qpA2R2 


CCACGACGGTATCTTGTGC 


68 


a 3 


qpA3Sl 


GCGAGATTCACGGTCCAATAA 


qpABRl 


GGCCATTTTGGTTTGTTTCG 


60 


a 4 


qpA4Sl 


GAGTATGGTGGCGTGCAAATG 


qpA4Rl 


GATATCCGGCCGCCAAAT 


60 


a 6 


qpA6Sl 


TGGAGAGACCTGTATCCAACGA 


qpA6Rl 


TGCTGTAGCGTGATGCCAAA 


64 


a7 


qpA7Sl 


CATGTATAATAGCGCTGACGAAGGT 


qpA7Rl 


CTGTTGACCACCACGTTGGTT 


63 


a 8 


qpASSl 


GAGGCACATCGACCAATCG 


qpASRl 


CGCTTAGATCAATGCCAACATC 


59 


a9 


qpA9Sl 


GTGCAACCCGTGCAGTACAG 


qpA9Rl 


TGCGTGTCATACGGCCAATA 


65 


a 10 


qpAlOSl 


GCACATGGTTCATAGCGAACTG 


qpAIORl 


GGTGTTCATATTCGCTCGGATT 


66 


P1 


qpBllSl 


CGCCGTCCAAACACAAGAT 


qpBllRl 


CTTGCAGTTGGGATGATGCA 


62 


P2 


qpB2Sl 


CCGTGAAGAGGAAAATACCG 


qpB2Rl 


GAACACGACGACTATCGCTG 


65 


rpl7 


qpRPL7F 


GCGCGCCGAGGCTTAT 


qpRPL7R 


CCGGAI 1 ILI 1 IGCAI 1 ILI IG 


81 


actin 


qpactinF 


AGCTCTATTCCAACCTTCCTTCT 


qpactinR 


TGTATGTAGTCTCGTGGATACCG 


62 



doi:l 0.1 371 /joumal.pone.0096669.t002 



Expression of nAChR subunits after neonicotinoid 
intoxication 

To study subunit expression levels after insecticide exposure, 
relative qPCR was performed on first larval stage LI exposed 
during 48 h with each neonicotinoid at LC50 or with DMSO 
(control condition). Total RNAs were extracted from intoxicated 
or control frrst larval stage (LI) using RNA Easy mini plant Kit 
and the same primers as described above. The results were 
normalized using the geometric mean of two reference genes, actin 
and rpl7 and validated using Normfmder software [22,23,24]. 
Expression levels were relative to control condition. 

Statistical analyses 

Statistical analyses were performed using GraphPad Prism 5 
(GraphPad Software Inc., La JoUa, CA). Data from binding 
experiments and insecticidal assays were analyzed by nonlinear 
regression analysis. A T-test (P<0.05; t-test with Welch's 
correction) was used for insecticidal assays and to compare Ki 
values. One-Way ANOVA (p<0.05) was used for qPCR 
experiments and binding assays. 

Results 

Binding properties of IMI, TMX and CLT on A.pisum native 
nAChRs 

Saturation binding experiments were carried out with both 
['^'^I]-a-Bgt and ['^H]-IMI on adult aphid membranes. Results are 
means of four experiments and the saturation binding parameters 
are summarized in Table 3. A saturation curve was first 
determined for [' ^'^L] -a-Bgt (Figure lA and IB) and revealed the 
presence of two binding sites: a high affinity (Kd = 
0.008±0.002 iiM and Bmax = 12.86±5.92 fmol/mg protein) 
and a low-affinity binding site (Kd = 1.135±0.213 nM and 
Bmax = 135.9±6.0 fmol/mg protein). The saturation data 
obtained for fH]-IMI (Figure 2 A and 2B) were also consistent 
with the presence of a high affinity (Kd = 0.16±0.04 nM and 
Bmax = 0.051 ±0.003 fmol/mg protein) and a low-affinity 
binding site (Kd = 4 1.7 ±5. 9 nM and Bmax = 
0.434±0.037 fmol/mg protein). For both ['ITI-IMI and ['^'^Fl-a- 



Bgt saturations curves, the presence of two binding sites was 
supported by the slope change in the Scatchard representation 
(Figure IB and 2B). In addition, we noticed that high affinity 
binding sites only represented 8.6±3.8% and 10.4±5.2% of total 
['^■^r|-a-Bgt and ['^IT|-IMI binding sites, respectively (Table 3). 
The difference between Bmax values for high- and low-affinity 
binding sites, for both ['^'iJ-a-Bgt and ['^H]-IMI, is consistent with 
the presence of these two binding sites on different nAChR 
populations in the aphid membrane preparation. In addition, the 
comparison of Bmax values highlighted a larger proportion of 
[ ' I]-a-Bgt binding sites compared to [' FT] -IMI (Table 3). Thus in 
the pea aphid A. pisum, ot-Bgt-sensitive nAChRs seem to represent 
a large majority of nAChR populations. 

In a second set of experiments, we studied the binding 
properties of IMI, TMX and CLT to the different nAChRs. 
Results are means of four experiments and the competition 
binding parameters are summarized in Table 3. For both 
radiotracers (['^H]-IMI and ['^"^FI-a-Bgt), we studied the low 
affinity binding sites, which represent the majority of iiACliR 
subtypes. Inhibition curves with IMI showed the presence of 20% 
of ['^'ri-a-Bgt residual binding in excess of IMI, suggesting that 
some of the a-Bgt binding sites were insensitive to IMI (Figure 3A). 
On the contrary, inhibition was complete using TMX (Figure 3B) 
and CLT (Figure 3C). Indeed, the inhibition constant (Ki) for ot- 
Bgt with a Kd value of 0.16 nM showed a better binding affinity 
for CLT and TMX (Ki = 0.18±0.05 pM and 1.53±0.65 ixM, 
respectively) compared to IMI (Ki = 14.61 ±1.13 jxM; One-Way 
ANOVA, p<0.05, table 3). With [''ITI-IMI, no specific binding 
inhibition was found using ot-Bgt, indicating that low- affinity ['^H] - 
IMI binding sites were insensitive to ot-Bgt (Figure 4A). Among the 
tested insecticides, homologous competition (Figure 4B) showed 
that IMI presents a Ki value of 38.14±6.88 nM which is 
consistent with the Kd value determined in the saturation 
experiment (41.7 nM, table 3). Interestingly, high concentrations 
of TMX and CLT were not able to completely displace [' H] -IMI 
from its binding sites, with maximal inhibition of 35% for TMX 
(Figure 4C) and 75% for CLT (Figure 4D). The apparent Ki 
values calculated from these binding curves were 1.05 ±0.07 |J,M 
for TMX and 127±42.5 iiM for CLT (Table 3). The residual 
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Figure 1. [^^^l]-a-Bungarotoxin specific binding on pea aphid. 

Saturation curves (A) and Scatchard plots (B) for ['^^l]-o(-Bungarotoxin 
(cz-Bgt) specific binding. Membranes were extracted from whole 
parthenogenetic adults of pea aphid Acyrthosiphon pisum LSR1. Results 
are means of four experiments. Error bars represent the SEM. 
doi:1 0.1 371 /journal.pone.0096669.g001 



Figure 2. [^H]-imidacloprid specific binding on pea aphid. 

Saturation curves (A) and Scatchard plots (B) for [^H]-imidacloprid 
specific binding. IVlembranes were extracted from whole parthenoge- 
netic adults of pea aphid Acyrthosiphon pisum LSR1. Results are means 
of four experiments. Error bars represent the SEM. 
doi:10.1371/journal.pone.0096669.g002 



binding could be explained by the inability of TMX and CLT to 
interact with all the nAChR subtypes recognized by IMI, or by an 
interaction of these ligands that were not strickly competitive with 
IMI. 



Toxicological effects of neonicotinoids on first-instar 
aphid larvae 

The toxicological effects of neonicotinoids upon A. pisum larvae 
have been previously studied [16]. Using the same method, we 
found that the three neonicotinoids TMX, CLT, and IMI have 



Table 3. [^^^l]-a-Bungarotoxin ([^^^IJ-a-Bgt) and pH]-imiclacloprid ([^H]-IMI) binding parameters determined on aphid membranes. 





Saturation binding 




['^'i]-a-Bgt 






[^Hl-IMI 








Kd (nM) 


Bmax (fmol/mg) 


% of high affinity 


Kd (nM) 


Bmax (fmol/mg) 


% of high affinity 


High affinity 


0.008±0.002 


12.86±5.92 


8.6±3.8 


0.16±0.04 


0.051 ±0.003 


10.4±5.2 


Low affinity 


1.135±0.213 


135 .9 ±6.0 




41.7±5.9 


0.434±0.037 




Competition binding 




['^=11- a-Bgt 






[^H]-IMI 








assay level (nM) 


IC50 (^M) 


Ki (nM) 


assay level (nM)IC50 (fiM) 


KI (nM) 


Imidacloprid 


0.6 


22.33 ±1.73 


14.6±1.13'' 


25 


0.061 ±0.011 


38.14±6.88= 


Thiamethoxam 


0.08 


1.636 ±0.701 


1.53±0.65'' 


25 


1.675±0.101 


1047.2±63.1'' 


Clothianidin 


0.08 


0.194±0.052 


0.18±0.05'' 


25 


0.203 ±0.068 


126.9±42.5= 


a-Bgt 








25 


n.d. 


n.d. 



Ki values were calculated according to Cheng and Prusoff formula considering Kd of low affinity binding sites, n.d: not determined. IC50: half maximal inhibitory 
concentration. Results are mean of four experiments and are represented ± SD. Ki values that are significantly different using One-Way ANOVA (p<0.05) are noted with 
different letters. 

doi:l 0.1 371/journal.pone.0096669.t003 
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A B C 




Figure 3. Neonicotinoids inhibition of [^^^l]-a-Bungarotoxin specific binding. Inhibition curves were determined on membranes of whole 
parthenogenetic adults of pea aphid Acyrthosiphon pisum for three neonicotinoids: A) imidacloprid (IMI), B) thiamethoxam (TMX) and C) clothianidin 
(CLT). Results are means of four experiments. Error bars represent the SEM. 
doi:1 0.1 371 /journal.pone.0096669.g003 



difTerent toxicities again.st A. pisum. All results are presented in 
table 4 as means of 6 to 8 experiments. We found that TMX was 
the most toxic (LC50 = 0.259 |xg/ml) and CLT was the least toxic 
(LC50 - 3.458 Hg/ml) after 24 h of exposure. The toxicity of IMI 
was intermediate with an LC-jo of 0.913 ng/ml. The LC50 values 
were significantly lower after 48 h of exposure. Interestingly IMI 
and TMX showed similar effect (LC50 = 0.038 and 0.034 Hg/ml, 
respectively) whereas CLT remained the least toxic (LC50 = 
0.118 |Xg/ml). We suggest that the potency of TMX could be 
associated to its double action: directly and after metabolization to 
CLT as previously proposed [25]. Indeed HPLC/MS analysis 
showed that TMX was metabolized to CLT (Figure 5). In TMX- 
treated aphids (using TMX at LC50 = 0.034 |j.g/ ml) we obtained 



final concentrations of 1.34 Hg/kg of TMX and 1.76 |J.g/kg of 
CLT, after 48 h of exposure. 

Is the expression of aphid nAChR subunits influenced by 
developmental stage or by exposure to neonicotinoids? 

Recently, using the full genome of A. pisum, Dale et al. 
highlighted the presence of 1 1 genes encoding putative nACliR 
subunits [15]. We confirmed the expression of these 11 genes in 
the pea aphid and studied the expression profile of these subunits 
according to the developmental stage and neonicotinoid exposure. 
First, qPCR experiments on the different developmental stages 
demonstrated that the expression of Apisumal, Apisuma2, 
Apisuma6, Apisuma8 and Apisum|32 was stable at the beginning 




log([CLT]) MM 



Figure 4. Neonicotinoid inhibition of [^H]-imidacioprid specific binding. Inhibition curves were determined on membranes of whole 
parthenogenetic adults of pea aphid Acyrthosiphon pisun) for (A) a-Bungarotoxin (cx-Bgt) and three neonicotinoids: (B) imidacloprid (IMI), (C) 
thiamethoxam (TMX) and (D) clothianidin (CLT). Results are means of four experiments. Error bars represent the SEM. 
doi:1 0.1 371 /journal.pone.0096669.g004 
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Figure 5. MRM chromatograms. Chromatograms of thiamethoxam (A) and its metabolite clothianidin (B) in 13,400 pea aphid larvae exposed to 
thiamethoxam at LC50 for 48 h. Intensity represents the peak area of the detected signal. 
doi:1 0.1 371/journal.pone.0096669.g005 



of aphid development and then was significantly reduced during 
adulthood (One Way Anova, p<0.05, n = 3 experiments in 
triplicate, figure 6). On the contrary ApisumaS expression 



increased with developmental stages although Apisuma?, Api- 
sumalO and ApisumPl transcript levels remained stable. The 
expression level of Apisumot4 and Apisum0(9 subunits showed 



Table 4. Neonicotinoid toxicity in the pea aphid on first-instar larvae for 24 h and 48 h of insecticide exposure using an artificial 
diet system. 







24 h 




48 h 






LC 50 (Mg/ml) 


CI 95% 


LC 50 (Mg/ml) 


CI 95% 


n 


Imidacloprid 


0.913 ' 


0.266 - 3.133 


0.038 


0.023 - 0.064 


2381 


Thiamethoxam 


0.259 


0.039 - 1.718 


0.034 


0.012 - 0.101 


2613 


Clothianidin 


3.458 " 


0.834 - 14.34 


0.118 " 


0.009 - 1 .62 


3016 



n = number of insects tested; CI = confidence interval; LC50 = Lethal concentration leading to 50% mortality. Results were corrected using Henderson-Tilton's formula. 
Values followed by different letters are significantly different {P<0.05; t-test with Welch's correction}. Toxicity curves were determined with 8 concentrations and 6 to 8 
replicates were made for each concentration. 
doi:l 0.1 371/journal.pone.0096669.t004 
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Figure 6. Expression level of nAChR mRNA subunits according to developmental stages of the pea aphid. Quantitative experiments 
were performed on vj/hole individuals in triplicate. Results are mean of three independent experiments. Relative expression ratio were calculated 
relative to first-instar nymphs and normalized with external reference gene luciferase. Statistical analysis (One-Way ANOVA) was carried out using 
GraphPad Prism 5 software. For each subunit, expression ratio statistically different according to larval stage are designated by different letters. 
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greater variability, with a lower expression level at the fourth larval 
stage. Thus, A. pisum subunits expression was regulated during 
developmental stages, suggesting that different nAChR subtypes 
could be expressed. 

Second, qPCR experiments were performed on surviving aphid 
larvae exposed to IMI, TMX and CLT after 48 h exposure. For 
this purpose, aphid larvae were intoxicated at the LC50 
determined in insecticidal assays. We found that IMI induced a 
strong variation of iiAChR subunits expression compared to 
control condition, with the exception of Apisumo(4, ApisumaS, 
Apisuma? and Apisuma9 (Figure 7, n = 4 to 7 experiments in 
triplicate). We also observed a significant increase of ApisumalO (+ 
218±40%), Apisumpi (+240±40%), Apisumal (+120±27%), 
Apisuma2 (+104±17%) and ApisumaS (+61±10%), respectively. 
On the contrary, a decrease was found with ApisumaS (— 34±4%) 
and ApisumP2 (— 40±4%), respectively (Figure 7A). Aphid larvae 
treated with TMX showed a significant decrease of Apisuma2 (— 
23±5%), Apisuma7 (-46±13%), Apisum|32 (-29±4%) and an 
increase for ApisumotlO (+90±31%) and Apisum|3l (+39±13%; 
figure 7B). Exposure to CLT led to a significant diminution of 
Apisuma4 (-49±4%), ApisumotS (-73±4%) and ApisumP2 (- 
48±3%) whereas we found a significant increase of ApisumalO (+ 
56± 12%; figure 7C). These data confirmed that the expression of 
aphid nAChR subunits was difFerentiaUy modified after exposure 
to various neonicotinoids. 

Discussion 

The pea aphid presents several pharmacological binding 
sites with different affinity for neonicotinoids 

In the present studies, saturation binding experiments demon- 
strated that fH]-IMI and ['^■"'Ij-Ol-Bgt labeled high- and low- 
affinity nAChR binding sites in A. pisum. Binding properties of ot- 
Bgt were not well documented in insects but two binding sites have 
also been found in the aphid M. persicae [13]. Interestingly, two 
specific [^FTj-IMI binding sites were likewise reported in the aphids 
M. persicae and A. craccivom [11,12]. In the pea aphid, the large 
difference between Bmax values for high and low affinity binding 
sites for both ['^FTj-IMI and ['^'Tj-ot-Bgt was in accordance with 
the presence of these sites on distinct nAChR subtypes and not at 
various subunit interfaces on the same receptor. In addition. 



competitive data showed that [■^H]-IMI low affinity binding site 
was insensitive to a-Bgt. We proposed that IMI could bind to ot- 
Bgt-insensitive nAChR subtypes, which was consistent with data 
obtained in M. persicae and A. craccivora [12]. Moreover, in the pea 
aphid, it seemed that the majority of binding sites was sensitive to 
a-Bgt, as previously demonstrated in D. melanogaster and M. persicae 
[11,12,13]. Competitive experiments also revealed that CLT 
bound to both ['^''I] -a-Bgt and [^H]-IMI binding sites. CLT was 
known to interact well with IMI-binding sites in the aphids M. 
persicae and A. craccivora [3,26]. Interestingly, only one study 
referred to competitive experiments between CLT and labeled a- 
Bgt. Zhang et al. demonstrated a weak inhibitor potency of CLT 
to [^H] -a-Bgt binding sites in D. melanogaster [27]. The apparent 
discrepancy between these results and ours could be attributed to 
species specificity. Furthermore, TMX, which was metabolized to 
CLT, showed a weak binding capacity for ['^H]-IMI binding sites 
and a better binding potency for ['^ 'Lj -a-Bgt binding sites in A. 
pisum. These results are consistent with previous studies describing 
a lack of TMX competition with [^H] -IMI in other aphid species 
such as M. persicae and A. craccivora [3,12]. Unfortunately, there was 
no data on TMX competitive binding to a-Bgt sites, despite that 
['H]-TMX could bind directiy in M. persicae and A. craccivora [28]. 
We propose that TMX binds to a-Bgt-binding sites in the pea 
aphid and that this mechanism could be present in other aphid 
species. 

The neonicotinoids IMI, TMX and CLT have different 
toxicological effects on A. pisum 

Acute toxicological assays demonstrated that TMX and IMI 
were more toxic than CLT. Similar data have been found with A. 
gossypii in which IMI was more toxic than both TMX and CLT 
[10]. Interestingly, in M. persicae CLT was found to be a more 
potent insecticide than IMI [8]. This discrepancy could be linked 
to variation in the intoxication method and/or susceptibility of 
aphid species. In other studies, the neonicotinoid susceptibility was 
evaluated using topical application and a dipping method for M. 
persicae and A. gossypii, respectively [8,10]. By contrast, we used an 
artificial diet protocol previously described by Sadeghi et al. [16]. 
The LC5Q at 48 h for IMI corresponded to the LC50 72 h in 
Sadeghi's study which indicated that the pea aphid strains could 
be differentially sensitive to neonicotinoids. Moreover, because 
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Figure 7. Expression levels of nAChR mRNA subunits after 
neonicotinoid exposure. Experiments were assessed on whole 
surviving! larvae exposed to neonicotinoids at LC50 for 48h. Aphids 
were intoxicated with imidacloprid (A) thiamethoxam (B) or clothianidin 
(C). Each qPCR experiment was performed in triplicate and results are 
represented as the mean of four to seven independent experiments 
after normalization with actin and ribosomal rp/7 gene. Error bars 
represent the SEM. Results are expressed in % of the expression level in 
control conditions (no insecticide, corresponding to 100%). Statistical 
analysis (t-test, a = 0.05) was carried out using Graphpad Prism 5 
software. 
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part of TMX was metabolized to CLT, we propose that the 
unusually high toxicity of TMX in A.pisum compared to other 
aphid species, was associated to its metabolite CLT [8,10]. This 
hypothesis has been previously demonstrated using the moth 
Spodopterafmgiperda and the cockroach Periplamta americana. Indeed, 
in S. Jmgiperda and P. americana, TMX was metabolized 24 h after 
treatment [25,29]. Nevertheless, in the present study, the 
proportion of metabolized TMX was different than previous 
studies [25,29]. Thus, the high insecticidal effect of TMX on the 
pea aphid could be due to its double action, by acting on uAChRs 
sensitive to ot-Bgt and IMI. 

Dale et al identified 1 1 putative genes encoding nAChR 
subunits in the pea aphid genome among which three were 
divergent (Apisum(32, Apisumot9 and ApisumalO) and did not 
belong to conserved subunit groups between insects species [15]. 
Using qPCR experiments on surviving larvae, we demonstrated 
that IMI, TMX and CLT significantly influenced iiAChR subunit 
expression. For all neonicotinoids tested, we found that Apisum 
alO was highly expressed after treatment. This subunit is an 



uncommon nAChR subunit lacking one cysteine in the Cys-loop 
and could be involved in distinct functional properties [15]. 
Previous studies performed with electric ray Torpedo demonstrated 
that a subunit lacking one cysteine in the Cys-loop could co- 
assemble to form functional receptors that are expressed at the 
membrane [30] . In Torpedo the lack of one cysteine also led to the 
loss of ot-Bgt binding sites [30] . As neonicotinoids bind to a-Bgt 
binding sites in the pea aphid, we propose that increased 
expression of ApisumalO subunit could likewise lead to increased 
expression of nAChR subtypes that are less sensitive to 
neonicotinoids. After neonicotinoid exposure, we also observed 
that Apisumfil was over-expressed after TMX and IMI exposure 
and Apisumfi2 under-expressed after treatment with the three 
insecticides. Thus, in the pea aphid we proposed that both 
Apisumfil and ApisumP2 could be differently involved in the 
regulation of neonicotinoid sensitivity. Indeed, studies performed 
on the brown planthopper Milapawata lugens and the aphid M. 
persicae showed that (31 was part of the IMI binding sites. Mutation 
of arginine to threonine at position 8 1 in this subunit induced an 
increase of resistance against neonicotinoids [31,32,33]. Moreover, 
recent studies demonstrated that nAChR subunit expression level 
was associated to neonicotinoid sensitivity [34,35] . For example, a 
decrease in AccPl and Accfi2 subunit expression in the Asiatic 
honey bee Apis cerana was described after IMI exposure [35]. Our 
results also demonstrated that neonicotinoid toxicity was associ- 
ated with specific iiAChR subunit regulation. For example, 
ApisumotlO and Apisum|32 were always up- or down-regulated, 
following treatment with any of the three insecticides. On the 
contrary, we found that Apisumot? and Apisuma4 were decreased 
only after exposure to TMX and CLT and the expression of 
Apisumot 1 and Apisuma3 was increased after IMI exposure. 
These results suggest that some subunits could be involved in 
specific insecticide action. Thus, we suggest that high and low 
affmity binding sites could involved several nAChR subtypes. 

Conclusions 

Previous studies conclude that nAChR subunits influence the 
pharmacological properties of nicotinic receptors and thus could 
modify the neonicotinoid sensitivity [36] . Our results demonstrat- 
ed that pea aphid iiAChR subunits were differentially expressed, 
first between developmental stages, as previously demonstrated in 
Drosophila and Apis cerana cerana [35,37,38], and also according to 
the neonicotinoid exposure. Neonicotinoid sensitivity could then 
be dependent on either physiological status and/ or environmental 
conditions in the pea aphid. Moreover, using toxicological and 
binding studies, we highhghted differences in neonicotinoid 
sensitivity in the pea aphid as compared to other aphid species 
and strains [8,10,16]. Thus, the insecticide strategies against aphid 
pests should be optimized for each particular species. In the pea 
aphid, the role of divergent subunits Apisuma9, ApisumalO and 
Apisumfi2 could be of particular interest to further understand the 
neonicotinoid mode of action. 
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